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Abstract 

We discuss the Anderson localization of electromagnetic fields in the 
fluctuating plasma induced by the gravitational density perturbation be- 
fore the recombination time of the Universe. Randomly distributed local- 
ized coherent electromagnetic fields emerge in the thermal equilibrium be- 
fore the recombination time. We argue that the localized coherent electric 
fields eventually produce cosmological magnetic fields after the decoupling 
time. 
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1 Introduction 



It has long been suspected that the presently observed galactic and cluster 
magnetic fields (lO^®5auss)[0 (and possibly the intercluster magnetic field « 
10^^ gauss, existence of which is under disputes |lll@||3|.) has a cosmological 
origin. There are two competing theories for the galactic magnetic field: the 
primeval scenario and the dynamo theory. The former explains the galactic 
magnetic field as a consequence of amplification of a seed magnetic field by 
protogalaxy collapse. On the other hand, in the more popular galactic dynamo 
theory g , the magnetic field has been exponentially amplified from a tiny 
seed magnetic field by differential galactic rotation and turbulence. One of the 
weak points of the primeval magnetic field scenario has been the lack of its 
viable explanation. 

Nevertheless, an alternative theory of density irregularities on the basis of 
the primeval magnetic field |^ is attractive, because it can explain the magni- 
tude of the density perturbation and the galactic magnetic field simultaneously. 
Recently a direct estimate of magnetic fields of young galaxies (z 2) has been 
performed by measuring the Faraday rotation of quasar spectra through the 
Lyman-a absorption systems to give 10^ gauss [|§. This large magnetic field 
of the galaxy in the early stage may revive interest in the primeval scenario. 
Further, even for the galactic dynamo theory we need a seed magnetic field, 
however small it is. 

In this paper we would like to point out that before the recombination time 
the cosmic plasma is strongly fluctuating at large scales due to the gravitational 
density perturbation. The plasma density fluctuation plays a role of randomly 
distributed scatterers for the electromagnetic fields and make it impossible for 
the electromagnetic fields to propagate if the wave length is longer than the 
elastic mean free path of electromagnetic wave by the scatterers. The localized 
modes are thermally excited and the electric field is more or less aligned in the 
localization region. The localized electric field induces a coherent current and 
then the magnetic field « 10^^® gauss at the time of decoupling of photons 
from matter. 

The localization phenomenon by random fields was first proposed by Ander- 
son ||l0|] to explain the behavior of the electric resistivity of metals in the presence 
of impurities and the theories were developed by himself and his collaborators 
[PI fl^. The localization phenomena are observed in various physical systems, 
e.g. in the propagation of light waves ||l^ and of sound waves jlj]. Probably 
the first paper which studied the possible localization of electromagnetic field in 
plasma is the one by Escande and Souillard JlSt . A preliminary experiment was 
performed for the longitudinal electron plasma wave by Dovcil ct al. (We 
shall adopt the localization theory in plasma ||l^] to the cosmological plasma in 
the expanding Universe.) 

Roughly speaking, the Anderson localization phenomenon occurs because of 
the destructive interference of the incident wave and the scattered waves from 
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the randomly distributed sources. It is known that the localization occurs if the 
strength of the randomness is larger than a critical value. For a sufficiently large 
random potential, only the modes of the wave length shorter than the mean free 
path can propagate, while the other modes are localized. 

We are going to argue that this localization phenomenon should occur in 
the matter dominated era in the standard cosmology. The resultant localized 
magnetic field may explain the origin of the galactic magnetic field and possibly 
a structure formation in the early stage |^ [Q . 

This paper is organized as follows. After briefly reviewing the electromag- 
netism in the cosmological plasma in §2, we will derive the basic equation for 
the Anderson localization of electromagnetic fields in the standard big bang 
model of the Universe in §3. In §4 and §5, we will estimate the mobility edge, 
the localization length and the number of localized states. We will describe a 
scenario which tells how the localized coherent electromagnetic fields emerge 
before the recombination time in §6 and give a rough estimate of the cosmolog- 
ical magnetic field at the decoupling time. The final section §7 is devoted to 
summary and discussions. 



2 The Plasma Equations in the Friedmann-Robertson- 
Walker Universe 

Here we will first review the plasma equations in the Fricdmann-Robertson- 
Walker background. 

Let us recall that the relativistic expression for the Maxwell equation is given 

by 

F''% = r, (1) 

where the electromagnetic field strength Fq,^ = daA/^ — d/sAa with being the 
vector potential. Here the current = enu" with u" being the four velocity 
of charged particles(electrons) with Uau" = —1 and n the density. The velocity 
u" satisfies the Lorentz force equation: 

mu"u\\a = Ff^aeu". (2) 

For simplicity, consider the spatially flat Friedmann- Robertson- Walker Universe 
with the metric 

ds^ ^aiXfi-dX^ +dx% (3) 

Here A is the conformal time defined by A = /* — . We assume that the charged 
particles are non-relativistic. So u" — (a^^,w*) with being small compared 
with unity. Then the Maxwell equation reads 

ri"^dpF,^ = ena^v,. (4) 
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with r]"'^ being the Lorentzian metric and the particle equation is approximately 
given by 

^~di^ eF,o, (5) 

where m is the electron mass and e is its charge. From the conservation of 
the current = 0, we see that na^ is approximately constant in our non- 
relativistic case before the recombination time. 



3 Random Potential by Fluctuating Plasma Den- 
sity 

Combining the equations and (|^), we can derive a general relativistic equa- 
tion for the vector potential A in terms of the conformal comoving coordinates 
(A,x) by 

Split the plasma density n into the homogeneous part n and the space dependent 
random part 6n; 

n — n + Sn. (7) 
Then the right hand side of the plasma equation becomes 

9 2 

, ,rec , ,rec 

[^^ + i^^]A, (8) 



where Wp"^^ = 

Here we have used the conservation law:n = rirec/o}^ before the recombina- 
tion time, with the scale factor a normalized at the recombination time, rirec is 
the electron density (therefore the baryon density) at the recombination time. 

The first term in the bracket is spatially constant but a decreasing function 
of the conformal time A. The second term represents the fluctuating random 
environment induced by the baryonic density perturbation — which is equal 

to ^ in our non-relativistic case. Knowing the standard theory of linear per- 
turbation in cosmology we see that ^ — — oc a in the matter dominated era, 
SO that the second term is independent of the conformal time in that era. We 

write the random potential as V^(x) = — which is a function of comoving 
coordinates x only. So far we have discussed the period after the equal time and 
before the recombination time and found that the random potential induced by 
the density perturbation is constant in time. As we shall see this is important 
for the stability of the Anderson localization. After the recombination time, the 
density of free electrons sharply decreases so that the random potential changes 
as rapidly as the other terms in the Maxwell equations so that the localization 
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breaks down after the recombination time. We will discuss this localization and 
de-localization phenomena in more detail in §6. 

For the matter dominated universe, the scale factor is given by a{t) = 
(^)2/3 = (3^)2 with A = / ^ = tVeUt^'^. Putting them all together 
we reduce the Maxwell equations to the following wave equation as 

-^ + S^ = [^'/'' + ^(")]^- 

Here = Qtl^JJ^^^ . 

We expand the vector potential A in terms of the complete orthonormal set 
{ipn} as 

A = ^a„/„(A)?A«(x)+c.c. (10) 

with a„ being arbitrary coefficients. Here ipn{^) is a normalized solution of the 
eigenvalue equation: 

[-A + nx)]V'«(x)=i?^^„(x). (11) 
The wave equation reduces to 

- 1^ + = 0, (12) 

with the normalization condition: i{fndxfn — fndxfn) — 1- 

Localization means that a number of modes {ipmU < no} are bound states 
exponentially falling off at infinity. The marginal state uq corresponds to the 
so-called mobility edge as we shall explain in the next section. 



4 Mobility Edge and Mean Free Path 

Throughout this paper we only consider the adiabatic baryonic density pertur- 
bation for simplicity, the spectrum of which is roughly like the picture below. To 
simplify our discussion we assume that the spectrum is sharply peaked around 
the Jeans length Ic at the equal time ^ and the distribution of perturbation is 
Gaussian. To understand the Anderson localization in the Universe before the 
recombination time, we have to evaluate the mobility edge of the wave number, 
which is roughly the inverse of the mean free path Imf of the electromagnetic 
wave by the random potential. Consider a Schrodinger equation for the spatial 
part of the electromagnetic field: 

[-A + F(x)]V' = E^i;, (13) 

•^We are fully aware that inclusion of the "scale invariant spectrum" for scales larger than 
Ic produces a hierarchy of structures. We would like to postpone the investigation in this 
direction to future works. 
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Figure 1: The spectrum of the adiabatic baryonic perturbation 



with X being the comoving coordinates. RecaU that the random potential 
y(x) = ujip'^'^'^^ does not depend on the conformal time in the matter dom- 
inated era. The characteristic length scale Ic of the density perturbation is 
given by the Jeans length at the equal timc^ which is roughly the same as 
the horizon scale at that time. After that time the length scale Ic develops as 
oc a. This gives Ic ~ lO^^cm at the recombination time, while the horizon is 
« 2 X W^cm. 

Then the mean free path is given by 

1^} = —. (14) 

where a is the elastic scattering cross section of the electromagnetic wave by 
a single potential and ric ~ is the density of "impurities". The elastic 

cross section is estimated roughly as « 47r/^ by the standard wave mechanics, 
because the "energy" is low and is roughly given hy E^l'^ ^ 1, while the potential 
is effectively gigantic V x ZJ^'^^ll « 5 x 10^^. (Here we have taken 

^ « 10~^ because the assumed baryonic perturbation can be comparable to 
the temperature fluctuation of CMB observed by COBE [0. However, the 
localization length and other results in the present work are insensitive to the 
magnitude of the density perturbation.) Therefore the localization should occur, 

^The equal time teq is the time when the matter and radiation energy density are equal. 
This value depends on the value of the present baryon density. We take "^^^g") = (tj:£c.\2/3 _ 

^(^eqj teg 

6 as a typical value for Qq = 0-2. 
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Figure 2: The characteristic length Ic, mean free path l^f and locahzation 
length ^. 

since the critical value of Vl"^ will be of order one. Then the mean free path is 
given by 



and the mobility edge is E* « ^j^, which is consistent with our low energy 
picture. 

5 Localization Length 

Fortunately in our present case, the randomness is huge so that an estimate of 



where Ic is the minimum length scale in which range the potential can be con- 
sidered as a constant (the same notation as the one by Escande and Souillard 
[p!5|). For the adiabatic baryonic perturbation, it will be reasonable to take 
the Jeans length at the equal time teq as the small scale cut-off Ic, which is 
~ IQ^^cm at the recombination time. In our present case, \V\/E'^ w 5 x 10'^^ so 
that ^ « lO^^cm at the recombination time. 

It is easy to estimate the number of localized states N by equating it with 
the number of states below the mobility edge E* — ^nl^K. Namely, 




(15) 




(16) 




(17) 
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We find that the density of the localized states is given by 

= yC/- (18) 

There will be around 2000 lumps of electromagnetic fields inside the horizon at 
the recombination time. We will not discuss its direct observational implica- 
tions in the present work, because there are uncertainties in the estimate of the 
mobility edge so that we should not take the predicted number of the lumps too 
seriously. 



6 Cosmological Magnetic Field 

The localized modes in the previous sections are excited in the thermal equi- 
librium after the equal time and before the recombination time. We estimate 
the magnitude of the thermally excited localized electric fields. The localized 
electric field E;oc is expanded as 

"0 

E(oc = - X] an/n(A)V'n(x) + c.c. (19) 
so that the thermal average of E^^^ is 

= E"''iQ7i^l/n|'IV'np. (20) 

Since Q/AT « 1, |/„|2 « Q/X and ~ we obtain 

K^^T/e, (21) 

in a localization region of the size ^. Here we have taken into ac;count the 
contribution from a single localized mode only, because the peaks of the other 
localized modes are randomly distributed an d wi ll be located outside of the 
region under consideration. The estimate of E^^^^ can also be derived by the 
equal partition law. Note that the propagating modes(n > no) contribute to 
the Stefan-Boltzmann law. 

At this stage wc realize that approximately T/ujp photons occupies a local- 
ized state. As far as T/ujp » 1, a coherent localized state of the size ^ is a 
good picture so that we can approximately describe it by a classical field: 

"0 

Ejoc = - X ^T/a;p/„(A)e„(x)V„(x) + c.c, (22) 

where e„(x) is the polarization vector, which is perpendicular to the gradient 
of the wave function ipn{^)- Therefore, in a localization domain of the size ^ 



7 



Figure 3: The diffusion of the localized electromagnetic field 



in the cosmic plasma, we have aligned electric and magnetic fields, which are 
perpendicular to each other and oscillate almost at the plasma frequency. 

One might worry about the damping of the magnetic structure of the size 
^ by photon diffusion. However, the diffusion does not occur to the localized 
modes. This is one of the virtues of the Anderson localization mechanism. 

Since the localization is a delicate interference phenomenon, the random 
potential T^(x) should vary slowly so that the motion of wave packet is not 
disturbed p^ . This condition is met before the recombination time, which we 
have discussed so far. After that time we have to take into account the recom- 
bination process in our equations for the electromagnetic fields. The number of 
free electrons sharply decreases, and the random potential changes in time. 

After recombination time the packets of the localized coherent electromag- 
netic fields begin to diffuse. We will see the time evolution of the packet by 
solving the Vlasov equation for the electron distribution and the Maxwell equa- 
tions for the electromagnetic fields. For our weak field we can linearize the 
Vlasov equation as 

^+v.fl + ^E.f^ = 0, (23) 
ot ox m ov 

where /i(t,x, v) is a small correction to the equilibrium Maxwell distribution 
function /o(v). We assume that the initial distribution function is the Maxwell 
distribution around the coherent electron velocity u — — ^E, 

/(t,X,v) =/o(v-u(x)) t<trec- (24) 

To solve the transport equation and the Maxwell equations, the standard method 
is to use the Laplace transformation, 

/>oo 

/i(^,v)= / dte'^'h{t,v). (25) 
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In the Laplace transformed space we can easily obtain an electric current density 
j = e / (Pvvfi, and then evaluate the Laplace inverse transform. Since the high 
frequency parts of the current density are averaged out, we evaluate only a part 
which originates from "the lu = kv pole" with k being the wave vector in the 
Fourier transformed space (see §34 in for the detailed calculation ). The 
result is ^ 

j(t, x) « -- / d3^/o(v)E(0, X - vt). (26) 

171 J 

In order to roughly estimate the total current /, we may simply suppose that the 
localized electromagnetic field has a Gaussian form E — (r/^^)^/^exp[— (x/^)^] 
at the recombination time. The induced magnetic B is then given by 



at the decoupling time with rg being the classical radius of electron. 

This linear picture holds until the v x B term becomes important, and after 
then the magnetohydrodynamics stage will set in. Since the localized region has 
a gigantic conductance, the magnetic field created at the decoupling era stays 
and will not dissipate to the Joule heat. 



7 Conclusion and Discussions 

In this paper we presented a new possible scenario which is intended to explain 
the origin of cosmological magnetic field. 

The baryonic density perturbation plays a role of random potential of the 
size lO^^cm to the electromagnetic fields and triggers the Anderson localization 
of the electric field before the recombination time. A large number of low 
energy photons occupy the localized states in the thermal equilibrium. The 
resultant coherent localized electric field is aligned in the localization region. 
The localized electric field then induces a coherent current and the magnetic 
field « 10~^^ gauss at the decoupling time. 

Roughly, we expect 2000 magnetic lumps of the size lO^^cm inside the hori- 
zon 2 X 10^'^ cm at that time. If the primeval magnetic field (10~^^ gauss) at 
the decoupling time is efficiently amplified by dynamo or other mechanism in 
the later stage, it may provide an explanation of the seed magnetic field in the 
Universe. However, it is beyond the scope of the present paper to work out 
the evolution of the magnetic fields in the long history of the Universe after 
the decoupling time. It is possible to argue that the resultant magnetic field 
traps sufficient number of protons so that the gravitational instability eventually 
creates a galaxy at the localization region. 

In the present work, we assumed that the large scale density perturbation 
is primordial and its magnitude is chosen to be around 4^ 10^^ so as to be 
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consistent with the COBE data on the anisotropy of CMB. In this case, the 
large scale structure is primordial and the galaxy scale one is its consequence by 
the Anderson localization. As remarked in the text, we only took into account 
a single scale Ic for the primordial perturbation for simplicity. The contribution 
from the larger scale to the localization phenomenon is interesting in the sense 
that it may account for the hierarchy of cosmological structures. 

It is well known that a density perturbation of a scale larger than the Jeans 
length I J = -jS= gravitationally collapses against the thermal pressure gradient 

VPG 



and therefore the perturbation grows. An analogous phenomenon exists for the 
gravitational collapse of localized magnetic fields. Correspondingly the magnetic 
Jeans length is obtained by replacing the sound velocity Cs in the expression of 
I J by the Alfven velocity Ca = B/^/p. 

^raJ = -7= (28) 



(For more satisfactory derivation, see Kim et al.0.) Just after recombination 
time, p sa lO^^'^.g/cm'^, we obtain 

Imj ~ — ^ X lO^^cm, (29) 
gauss 

which is comfortably less than the previously estimated localization length 
lO'^^cm, if B is less than the present galactic value 10^^ . gauss. Therefore the 
gravity exceeds the magnetic pressure so that the magnetic lumps will not blow 
up after the decoupling time. 
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Figure captions 



• Figure 1. The spectrum of the adiabatic baryonic perturbation. The dashed 
and soUd Unes represent the initial scale invariant spectrum and the spectrum 
at the recombination time, respectively. Ic is the Jeans length at the equal time. 

• Figure 2. The diagram shows the time evolution of some characteristic 
lengths, the Hubble horizon, the Jeans length at equal time Ic, the mean free 
path l^f, and the localization length ^. 

• Figure 3. The diffusion of the localized electromagnetic field. The upper 
figure represents the localized field before the recombination time, while the 
lower one represents the diffusion process after that time. 



12 



